Adhesion energy controls lipid binding-mediated endocytosis.

Several bacterial toxins and viruses can deform membranes through multivalent binding to lipids for clathrin-independent endocytosis. However, it remains unclear, how membrane deformation and endocytic internalization are mechanistically linked. Here we show that many lipid-binding virions induce membrane deformation and clathrin-independent endocytosis, suggesting a common mechanism based on multivalent lipid binding by globular particles. We create a synthetic cellular system consisting of a lipid-anchored receptor in the form of GPI-anchored anti-GFP nanobodies and a multivalent globular binder exposing 180 regularly-spaced GFP molecules on its surface. We show that these globular, 40 nm diameter, particles bind to cells expressing the receptor, deform the plasma membrane upon adhesion and become endocytosed in a clathrin-independent manner. We explore the role of the membrane adhesion energy in endocytosis by using receptors with affinities varying over 7 orders of magnitude. Using this system, we find that once a threshold in adhesion energy is overcome to allow for membrane deformation, endocytosis occurs reliably. Multivalent, binding-induced membrane deformation by globular binders is thus sufficient for internalization to occur and we suggest it is the common, purely biophysical mechanism for lipid-binding mediated endocytosis of toxins and pathogens.

Unravelling the mechanism of neurotensin recognition by neurotensin receptor 1.

The conformational ensembles of G protein-coupled receptors (GPCRs) include inactive and active states. Spectroscopy techniques, including NMR, show that agonists, antagonists and other ligands shift the ensemble toward specific states depending on the pharmacological efficacy of the ligand. How receptors recognize ligands and the kinetic mechanism underlying this population shift is poorly understood. Here, we investigate the kinetic mechanism of neurotensin recognition by neurotensin receptor 1 (NTS1) using 19F-NMR, hydrogen-deuterium exchange mass spectrometry and stopped-flow fluorescence spectroscopy. Our results indicate slow-exchanging conformational heterogeneity on the extracellular surface of ligand-bound NTS1. Numerical analysis of the kinetic data of neurotensin binding to NTS1 shows that ligand recognition follows an induced-fit mechanism, in which conformational changes occur after neurotensin binding. This approach is applicable to other GPCRs to provide insight into the kinetic regulation of ligand recognition by GPCRs.

Adhesion energy controls lipid binding-mediated endocytosis.

Several bacterial toxins and viruses can deform membranes through multivalent binding to lipids for clathrin-independent endocytosis. However, it remains unclear, how membrane deformation and endocytic internalization are mechanistically linked. Here we show that many lipid-binding virions induce membrane deformation and clathrin-independent endocytosis, suggesting a common mechanism based on multivalent lipid binding by globular particles. We create a synthetic cellular system consisting of a lipid-anchored receptor in the form of GPI-anchored anti-GFP nanobodies and a multivalent globular binder exposing 180 regularly-spaced GFP molecules on its surface. We show that these globular, 40 nm diameter, particles bind to cells expressing the receptor, deform the plasma membrane upon adhesion and become endocytosed in a clathrin-independent manner. We explore the role of the membrane adhesion energy in endocytosis by using receptors with affinities varying over 7 orders of magnitude. Using this system, we find that once a threshold in adhesion energy is overcome to allow for membrane deformation, endocytosis occurs reliably. Multivalent, binding-induced membrane deformation by globular binders is thus sufficient for internalization to occur and we suggest it is the common, purely biophysical mechanism for lipid-binding mediated endocytosis of toxins and pathogens.

Measuring pico-Newton Forces with Lipid Anchors as Force Sensors in Molecular Dynamics Simulations.

Binding forces between biomolecules are ubiquitous in nature but sometimes as weak as a few pico-Newtons (pN). In many cases, the binding partners are attached to biomembranes with the help of a lipid anchor. One important example are glycolipids that promote membrane adhesion through weak carbohydrate-carbohydrate binding between adjacent membranes. Here, we use molecular dynamics (MD) simulations to quantify the forces generated by bonds involving membrane-anchored molecules. We introduce a method in which the protrusion of the lipid anchors from the membrane acts as the force sensor. Our results with two different glycolipids reveal binding forces of up to 20 pN and corroborate the recent notion that carbohydrate-carbohydrate interactions are generic rather than specific.

Molecular Dynamics Simulations of Immune Receptors and Ligands.

Molecular dynamics simulations of immune receptor and ligand proteins in their native membrane environment allow to determine the orientational and structural variability of the proteins and protein complexes. The simulations complement the static, "membrane-free" structural information obtained from cryo-EM structures of transmembrane proteins in detergent micelles or from crystal structures of extracellular protein domains. Here we describe how to set up and perform simulations of transmembrane receptors, ligands, and receptor-ligand complexes.

Mechanical forces impair antigen discrimination by reducing differences in T-cell receptor/peptide-MHC off-rates.

T cells use their T-cell receptors (TCRs) to discriminate between lower-affinity self and higher-affinity foreign peptide major-histocompatibility-complexes (pMHCs) based on the TCR/pMHC off-rate. It is now appreciated that T cells generate mechanical forces during this process but how force impacts the TCR/pMHC off-rate remains debated. Here, we measured the effect of mechanical force on the off-rate of multiple TCR/pMHC interactions. Unexpectedly, we found that lower-affinity TCR/pMHCs with faster solution off-rates were more resistant to mechanical force (weak slip or catch bonds) than higher-affinity interactions (strong slip bonds). This was confirmed by molecular dynamics simulations. Consistent with these findings, we show that the best-characterized catch bond, involving the OT-I TCR, has a low affinity and an exceptionally fast solution off-rate. Our findings imply that reducing forces on the TCR/pMHC interaction improves antigen discrimination, and we suggest a role for the adhesion receptors CD2 and LFA-1 in force-shielding the TCR/pMHC interaction.

A litmus test for classifying recognition mechanisms of transiently binding proteins.

Partner recognition in protein binding is critical for all biological functions, and yet, delineating its mechanism is challenging, especially when recognition happens within microseconds. We present a theoretical and experimental framework based on straight-forward nuclear magnetic resonance relaxation dispersion measurements to investigate protein binding mechanisms on sub-millisecond timescales, which are beyond the reach of standard rapid-mixing experiments. This framework predicts that conformational selection prevails on ubiquitin's paradigmatic interaction with an SH3 (Src-homology 3) domain. By contrast, the SH3 domain recognizes ubiquitin in a two-state binding process. Subsequent molecular dynamics simulations and Markov state modeling reveal that the ubiquitin conformation selected for binding exhibits a characteristically extended C-terminus. Our framework is robust and expandable for implementation in other binding scenarios with the potential to show that conformational selection might be the design principle of the hubs in protein interaction networks.

Interplay of Trans- and Cis-Interactions of Glycolipids in Membrane Adhesion.

Glycolipids mediate stable membrane adhesion of potential biological relevance. In this article, we investigate the trans- and cis-interactions of glycolipids in molecular dynamics simulations and relate these interactions to the glycolipid-induced average separations of membranes obtained from neutron scattering experiments. We find that the cis-interactions between glycolipids in the same membrane leaflet tend to strengthen the trans-interactions between glycolipids in apposing leaflets. The trans-interactions of the glycolipids in our simulations require local membrane separations that are significantly smaller than the average membrane separations in the neutron scattering experiments, which indicates an important role of membrane shape fluctuations in glycolipid trans-binding. Simulations at the experimentally measured average membrane separations provide a molecular picture of the interplay between glycolipid attraction and steric repulsion of the fluctuating membranes probed in the experiments.

Structural variability and concerted motions of the T cell receptor - CD3 complex.

We investigate the structural and orientational variability of the membrane-embedded T cell receptor (TCR) - CD3 complex in extensive atomistic molecular dynamics simulations based on the recent cryo-EM structure determined by Dong et al., 2019. We find that the TCR extracellular (EC) domain is highly variable in its orientation by attaining tilt angles relative to the membrane normal that range from 15° to 55°. The tilt angle of the TCR EC domain is both coupled to a rotation of the domain and to characteristic changes throughout the TCR - CD3 complex, in particular in the EC interactions of the Cß FG loop of the TCR, as well as in the orientation of transmembrane helices. The concerted motions of the membrane-embedded TCR - CD3 complex revealed in our simulations provide atomistic insights on conformational changes of the complex in response to tilt-inducing forces on antigen-bound TCRs.

Cooperative Stabilization of Close-Contact Zones Leads to Sensitivity and Selectivity in T-Cell Recognition.

T cells are sensitive to 1 to 10 foreign-peptide-MHC complexes among a vast majority of self-peptide-MHC complexes, and discriminate selectively between peptide-MHC complexes that differ not much in their binding affinity to T-cell receptors (TCRs). Quantitative models that aim to explain this sensitivity and selectivity largely focus on single TCR/peptide-MHC complexes, but T cell adhesion involves a multitude of different complexes. In this article, we demonstrate in a three-dimensional computational model of T-cell adhesion that the cooperative stabilization of close-contact zones is sensitive to one to three foreign-peptide-MHC complexes and occurs at a rather sharp threshold affinity of these complexes, which implies selectivity. In these close-contact zones with lateral extensions of hundred to several hundred nanometers, few TCR/foreign-peptide-MHC complexes and many TCR/self-peptide-MHC complexes are segregated from LFA-1/ICAM-1 complexes that form at larger membrane separations. Previous high-resolution microscopy experiments indicate that the sensitivity and selectivity in the formation of closed-contact zones reported here are relevant for T-cell recognition, because the stabilization of close-contact zones by foreign, agonist peptide-MHC complexes precedes T-cell signaling and activation in the experiments.

Membrane morphologies induced by mixtures of arc-shaped particles with opposite curvature.

Biological membranes are shaped by various proteins that either generate inward or outward membrane curvature. In this article, we investigate the membrane morphologies induced by mixtures of arc-shaped particles with coarse-grained modeling and simulations. The particles bind to the membranes either with their inward, concave side or their outward, convex side and, thus, generate membrane curvature of opposite sign. We find that small fractions of convex-binding particles can stabilize three-way junctions of membrane tubules, as suggested for the protein lunapark in the endoplasmic reticulum of cells. For comparable fractions of concave-binding and convex-binding particles, we observe lines of particles of the same type, and diverse membrane morphologies with grooves and bulges induced by these particle lines. The alignment and segregation of the particles is driven by indirect, membrane-mediated interactions.

Weak carbohydrate-carbohydrate interactions in membrane adhesion are fuzzy and generic.

Carbohydrates such as the trisaccharide motif LeX are key constituents of cell surfaces. Despite intense research, the interactions between carbohydrates of apposing cells or membranes are not well understood. In this article, we investigate carbohydrate-carbohydrate interactions in membrane adhesion as well as in solution with extensive atomistic molecular dynamics simulations that exceed the simulation times of previous studies by orders of magnitude. For LeX, we obtain association constants of soluble carbohydrates, adhesion energies of lipid-anchored carbohydrates, and maximally sustained forces of carbohydrate complexes in membrane adhesion that are in good agreement with experimental results in the literature. Our simulations thus appear to provide a realistic, detailed picture of LeX-LeX interactions in solution and during membrane adhesion. In this picture, the LeX-LeX interactions are fuzzy, i.e. LeX pairs interact in a large variety of short-lived, bound conformations. For the synthetic tetrasaccharide Lac 2, which is composed of two lactose units, we observe similarly fuzzy interactions and obtain association constants of both soluble and lipid-anchored variants that are comparable to the corresponding association constants of LeX. The fuzzy, weak carbohydrate-carbohydrate interactions quantified in our simulations thus appear to be a generic feature of small, neutral carbohydrates such as LeX and Lac 2.

Accessory mutations balance the marginal stability of the HIV-1 protease in drug resistance.

The HIV-1 protease is a major target of inhibitor drugs in AIDS therapies. The therapies are impaired by mutations of the HIV-1 protease that can lead to resistance to protease inhibitors. These mutations are classified into major mutations, which usually occur first and clearly reduce the susceptibility to protease inhibitors, and minor, accessory mutations that occur later and individually do not substantially affect the susceptibility to inhibitors. Major mutations are predominantly located in the active site of the HIV-1 protease and can directly interfere with inhibitor binding. Minor mutations, in contrast, are typically located distal to the active site. A central question is how these distal mutations contribute to resistance development. In this article, we present a systematic computational investigation of stability changes caused by major and minor mutations of the HIV-1 protease. As most small single-domain proteins, the HIV-1 protease is only marginally stable. Mutations that destabilize the folded, active state of the protease therefore can shift the conformational equilibrium towards the unfolded, inactive state. We find that the most frequent major mutations destabilize the HIV-1 protease, whereas roughly half of the frequent minor mutations are stabilizing. An analysis of protease sequences from patients in treatment indicates that the stabilizing minor mutations are frequently correlated with destabilizing major mutations, and that highly resistant HIV-1 proteases exhibit significant fractions of stabilizing mutations. Our results thus indicate a central role of minor mutations in balancing the marginal stability of the protease against the destabilization induced by the most frequent major mutations.

On the relationship between docking scores and protein conformational changes in HIV-1 protease.

We present a detailed investigation of the effect of the crystallographic structure of the HIV-1 protease (PR) on the binding energy of different classes of inhibitors obtained from docking simulations. The crystal structures of 222 HIV-1 proteases (in wild-type and mutant forms) and 202 inhibitors were downloaded from appropriate databases. A cross-docking approach (docking of all 202 inhibitors to all 222 PR structures) using Autodock Vina was implemented. The protease structures were clustered using a Kohonen self-organization map analysis of the data matrix of docking energies. The obtained clusters of PRs were correlated with the x-y-z coordinates of the PR structures to identify structural segments underlying this clustering. The PR structures were clustered into 4 classes. One of these classes exhibits rather strong binding with almost all inhibitors, while another class exhibits rather weak binding. The remaining two classes are intermediate in binding strength. The selectivity ratio indices for the carbon-alpha atoms of the PR structures indicate that conformational motion of residues outside the binding pocket contributes significantly to the discrimination of the 4 classes.

Membrane Morphologies Induced by Arc-Shaped Scaffolds Are Determined by Arc Angle and Coverage.

The intricate shapes of biological membranes such as tubules and membrane stacks are induced by proteins. In this article, we systematically investigate the membrane shapes induced by arc-shaped scaffolds such as proteins and protein complexes with coarse-grained modeling and simulations. We find that arc-shaped scaffolds induce membrane tubules at membrane coverages larger than a threshold of ∼40%, irrespective of their arc angle. The membrane morphologies at intermediate coverages below this tubulation threshold, in contrast, strongly depend on the arc angle. Scaffolds with arc angles of about 60°, akin to N-BAR domains, do not change the membrane shape at coverages below the tubulation threshold, whereas scaffolds with arc angles larger than about 120° induce double-membrane stacks at intermediate coverages. The scaffolds stabilize the curved membrane edges that connect the membrane stacks, as suggested for complexes of reticulon proteins. Our results provide general insights on the determinants of membrane shaping by arc-shaped scaffolds.

Membrane fluctuations and acidosis regulate cooperative binding of 'marker of self' protein CD47 with the macrophage checkpoint receptor SIRPα.

Cell-cell interactions that result from membrane proteins binding weakly in trans can cause accumulations in cis that suggest cooperativity and thereby an acute sensitivity to environmental factors. The ubiquitous 'marker of self' protein CD47 binds weakly to SIRPα on macrophages, which leads to accumulation of SIRPα (also known as SHPS-1, CD172A and SIRPA) at phagocytic synapses and ultimately to inhibition of engulfment of 'self' cells - including cancer cells. We reconstituted this macrophage checkpoint with GFP-tagged CD47 on giant vesicles generated from plasma membranes and then imaged vesicles adhering to SIRPα immobilized on a surface. CD47 diffusion is impeded near the surface, and the binding-unbinding events reveal cooperative interactions as a concentration-dependent two-dimensional affinity. Membrane fluctuations out-of-plane link cooperativity to membrane flexibility with suppressed fluctuations in the vicinity of bound complexes. Slight acidity (pH 6) stiffens membranes, diminishes cooperative interactions and also reduces 'self' signaling of cancer cells in phagocytosis. Sensitivity of cell-cell interactions to microenvironmental factors - such as the acidity of tumors and other diseased or inflamed sites - can thus arise from the collective cooperative properties of flexible membranes.This article has an associated First Person interview with the first author of the paper.

Membrane-Mediated Cooperativity of Proteins.

Besides direct protein-protein interactions, indirect interactions mediated by membranes play an important role for the assembly and cooperative function of proteins in membrane shaping and adhesion. The intricate shapes of biological membranes are generated by proteins that locally induce membrane curvature. Indirect curvature-mediated interactions between these proteins arise because the proteins jointly affect the bending energy of the membranes. These curvature-mediated interactions are attractive for crescent-shaped proteins and are a driving force in the assembly of the proteins during membrane tubulation. Membrane adhesion results from the binding of receptor and ligand proteins that are anchored in the apposing membranes. The binding of these proteins strongly depends on nanoscale shape fluctuations of the membranes, leading to a fluctuation-mediated binding cooperativity. A length mismatch between receptor-ligand complexes in membrane adhesion zones causes repulsive curvature-mediated interactions that are a driving force for the length-based segregation of proteins during membrane adhesion.

Author Correction: Protein-peptide association kinetics beyond the seconds timescale from atomistic simulations.

In the original version of this Article, the Acknowledgement section omitted financial support from the Deutsche Forschungsgemeinschaft grant SFB 958/A4. This error has now been corrected in both the PDF and HTML versions of the Article.

Identifying Conformational-Selection and Induced-Fit Aspects in the Binding-Induced Folding of PMI from Markov State Modeling of Atomistic Simulations.

Unstructured proteins and peptides typically fold during binding to ligand proteins. A challenging problem is to identify the mechanism and kinetics of these binding-induced folding processes in experiments and atomistic simulations. In this Article, we present a detailed picture for the folding of the inhibitor peptide PMI into a helix during binding to the oncoprotein fragment 25-109Mdm2 obtained from atomistic, explicit-water simulations and Markov state modeling. We find that binding-induced folding of PMI is highly parallel and can occur along a multitude of pathways. Some pathways are induced-fit-like with binding occurring prior to PMI helix formation, while other pathways are conformational-selection-like with binding after helix formation. On the majority of pathways, however, binding is intricately coupled to folding, without clear temporal ordering. A central feature of these pathways is PMI motion on the Mdm2 surface, along the binding groove of Mdm2 or over the rim of this groove. The native binding groove of Mdm2 thus appears as an asymmetric funnel for PMI binding. Overall, binding-induced folding of PMI does not fit into the classical picture of induced fit or conformational selection that implies a clear temporal ordering of binding and folding events. We argue that this holds in general for binding-induced folding processes because binding and folding events in these processes likely occur on similar time scales and do exhibit the time-scale separation required for temporal ordering.